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Summary. In the present paper different types of scramjet (supersonic combus¬ 
tion ramjet) combustors are investigated. Thereby the main difference between the 
combustors is the way of injecting the fuel into the combustion chamber. The first 
investigated concept of fuel injection is the injection by strut injectors. Here the 
injection of fuel is realized by a lobed strut that is located in the middle of the 
combustion chamber. The second concept for fuel supply is the wall injection of 
hydrogen. Here the fuel is injected by several holes in the wall of the combustor. 
Both concepts of fuel injection have different advantages and disadvantages which 
are explained in detail. Although different performance parameters for both scram¬ 
jet combustors are introduced this paper will not compare the different techniques 
among each other. Because of the high Reynolds numbers in scramjet combustors, 
the need to resolve the boundary layers and the necessity of detailed chemistry, the 
simulation of scramjets is extremely CPU time demanding. 


1 Introduction 

Due to the high velocities in a scramjet combustor the residence time of air and 
fuel in the combustion chamber is extremely short. Additionaly the mixing 
rates at high flow Mach numbers are inherently low. Hence techniques are 
required to achieve a rapid and efficient mixing of fuel and air and thus a 
stable combustion and a complete burnout of the fuel. There are mainly two 
concepts of fuel injection: 1. wall injectors [1, 2, 3] and 2. strut injectors 
[4, 5, 6, 7]. In case of a wall injection the fuel is injected through the wall into 
the air flow. Wall injectors are easy to manufacture and easy to cool, have a 
good near field mixing and cause no pressure loss when they are turned off. 
On the other hand in real size combustors there might be problems of the 
penetration depth of the fuel into the air flow and a strong blockage might 
occur in the combustion chamber. Strut injectors, where the fuel is injected 
through a strut directly into the core of the air flow, cannot be removed from 
the flow and hence cause pressure losses even when they are switched off. 
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Additionally they have to be cooled and usually need mixing enhancement 
techniques. The advantages of strut injectors are the absence of strong shock 
waves due to a blockage caused by the fuel jet. In the present paper both 
concepts are investigated numerically. 


2 Governing Equations and Numerical Scheme 


The scientific inhouse code TASCOM3D (Turbulent All Speed Combustion 
Multigrid) describes reacting flows by solving the full compressible Navier- 
Stokes, species and turbulence transport equations. Additionally an assumed 
PDF (probability density function) approach is used to take turbulence chem¬ 
istry interaction into consideration. Therefore two additional equations (for 
the variance of the temperature and the variance of the sum of species mass 
fractions) have to be solved. Thus the described set of averaged equations in 
three-dimensional conservative form is given by 
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The variables in the conservative variable vector Q are the density p (aver¬ 
aged), the velocity components (Favre averaged) u, v and w, the total specific 
energy E, the turbulence variables q = \/k and lo = e/k (where k is the kinetic 
energy and e the dissipation rate of k), the variance of the temperature <jt 
and the variance of the sum of the species mass fractions ay and finally the 
species mass fractions (* = 1,2,..., AT fe — 1). Thereby N k describes the total 
number of species that are used for the description of the gas composition. 
The vectors F, G and H specify the inviscid fluxes in x-, y- and z- direction, 
F„, G„ and H„ the viscous fluxes, respectively. The source vector S in Eq. 
(1) results from turbulence and chemistry and is given by 

S=[0,0,0,0,0,S q ,S u ,S <TT ,S <TY ,S Yi ] T , i = l,2,...,N k -l, (3) 

where S q and S u are the averaged source terms of the turbulence variables, 
S aT and S,y Y the source terms of the variance variables (<tt and ay) and 
Sy i the source terms of the species mass fractions. For turbulence closure a 
two-equation low-Reynolds-number q-u turbulence-model is applied [8]. The 
momentary chemical production rate of species i in Eq. (1) is defined by 
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where kf r and /q v are the forward and backward rate constants of reaction 
r (defined by the Arrhenius function), the molecular weight of a species Mi, 
the species concentration Ci = pY) /Mj and the stoichiometric coefficients z/, ; r 
and v i r of species i in reaction r. The averaged chemical production rate for a 
species i due to the use of the assumed PDF approach and the assumption of 
statistical independence of temperature, gas composition and density is given 

by 

Sy, = J Si ■■ , c Nk ^j P (f,c i, • • •, cjv fe ) dTdci ■ ■ ■ dc Nk , (5) 

where 

P {t , cp • • •, c Nk ) = P T (f) Py (n, • ■ ■, Y Nk ) S(p-p) . (6) 


Pt defines the temperature PDF desribed by a Gaussian distribution 
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that is clipped at lower and upper temperature limits due to the limitations of 
the Arrhenius equation [9]. The PDF of the gas composition Py is described 
by the multi-variate /3-pdf proposed by Girimaji [10] 
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5 in Eq. (6) and (8) denotes the (5-function. The unsteady set of equations (1) is 
solved using an implicit Lower-Upper-Symmetric Gauss-Seidel (LU-SGS) [11, 
12, 13, 14] finite-volume algorithm, where the finite-rate chemistry is treated 
fully coupled with the fluid motion. More details concerning TASCOM3D may 
be found in [13, 14, 15, 16, 17, 18] 


3 Investigations of Lobed Strut Injectors 

In this section a scramjet combustion chamber with a strut injector is investi¬ 
gated. Main focus is the effect of different amounts of circulation produced by 
the lobed strut on the mixing and hence on the combustion process. Thereby 
hydrogen is injected in axial flow direction by a lobed strut injector (see Fig. 1) 
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which creates streamwise vortices to enhance the mixing (compared to planar 
struts) [19]. In order to produce different amounts of circulation the length of 
the lobed structure of the strut is varied. Fig. 1 shows the struts investigated 
which all have 60 mm in width and 8 mm in height, while the length of the 
lobed parts differ. The area for hydrogen injection is the same in all three 
cases. 


strut I 




Fig. 1 . Investigated struts with different lengths and amount of vorticity production 



Fig. 2. Sketch of the combustion chamber with a lobed strut injector 


In case of strut I the length is 80 mm. For strut II the lobed part is short¬ 
ened by a factor of 2 (60 mm total length) and for strut III by a factor of 
4 in comparison to strut I (50 mm total length). A reduction of the strut 
length increases the angles for changes in flow direction and thus the strength 
of the expansion fans at the middle of the strut and the pressure ratio in 
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cross stream direction. Thus the shorter struts create stronger vortices. Fig. 
2 demonstrates a sketch of the combustion chamber with a lobed strut in¬ 
jector. The dimensions of the combustor are 536 mm x 60 mm x 38 mm. 
Thereby the computional domain includes only a section of the combustion 
chamber using periodic and symmetry boundary comditions. Furthermore the 
simulations have been performed in two steps: a two-dimensional simulation 
up to the middle of the strut and three-dimensional simulations of the re¬ 
maining combustion chamber. The inflow conditions for the simulations are 
summarized in Tab. 1. More details concerning the model scramjet combus¬ 
tor geometry and the lobed strut injectors may be found in [20]. The grid of 
the computational domain has about 800.000 volumes (220 x 70 x 54) and 
is strongly refined at all near wall regions as well as in the main combustion 
zone. Although the quantity of the volumes is moderate the computational 
costs are high due to the additional transport equations for the variances 
and especially the species as well as the numerical stiffness caused by com¬ 
bustion. Fig. 3 shows representatively for all configurations the H 2 O, H 2 and 
OH distributions of the model scramjet combustor using strut I, respectively. 
The hydrogen distribution nicely shows the production of vortices due to the 
lobed strut injector right after the injection. The H 2 O and OH distributions 
represent the combustion process and show that in all cases (strut I - strut 
III) lifted flames are obtained. As observed in previous investigations [19], 
ignition takes place near the point where the diverging channel part begins. 
The positions for ignition vary between 85 mm (strut I), 105 mm (strut II) 
and 110 mm (strut III), respectively. Further downstream the flame becomes 
wider and the separated structures grow together till the flame has a W-shape 
at the outlet of the combustion chamber. In order to describe the effect of the 
struts on the combustion process calculated H 2 O distributions are plotted in 
Fig. 4 for different y-z slices at channel positions ranging from 100-400 mm. 
100 mm downstream of the hydrogen injection the shortest strut has a very 
small amount of H 2 0, because ignition has just occured. In case of the longest 


Table 1. Inflow conditions for the simulation of the combustion chamber with the 
lobed strut injector 



air 

strut 

pressure p [Pa] 

211000 

50000 

temperature T [K] 

1300 

196 

velocity u [™] 

723 

2281 

Mach-number Ma [-] 

1 

2 

mass fraction Ytt 
h 2 

0 

1 

mass fraction Yo 
u 2 

0.23 

0 

mass fraction Ym 

.. z 

0.77 

0 

m [f] 

628 

4.18 
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Fig. 3. Calculated H 2 , H 2 O and OH distributions (from top to bottom) of the model 
scramjet combustor using strut I, respectively 


strut the combustion process is already in progress and the amount of H 2 O is 
much higher compared to the other struts. The main combustion zone in all 
three cases forms two circles which are turning to a more elliptical shape with 
an increased strut length. With increasing distance from the fuel injector the 
flame spreads in vertical direction and the circle-like shapes grow together. In 
case of the shortest strut the flame shows the most homogeneous H 2 O distri¬ 
bution, whereas the differences are quite small at the end of the combustion 
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H,0= 0.02 0.06 0.1 0.14 0.18 0.22 0.26 0.3 


Fig. 4. Calculated H 2 O molar fraction in y-z cross sections at x = 100, 200, 300, 
and 400 mm for the struts I to III (from top to bottom) 


chamber. Finally in Fig. 5 a comparison of the three struts by means of per¬ 
formance parameters (hydrogen mass flux, total pressure) is given. Strut I has 
the best mixing in the near field of the injector which causes an early ignition 
and burning of the hydrogen (indicated by the rapid decrease of the hydrogen 
mass flux at the normalized channel length of 0.15). However, with increas¬ 
ing distance from the strut, vorticity seems to be too weak to transport air 
from the upper and lower channel walls into the main combustion zone. This 
results in a higher amount of unburned hydrogen at the end of the combustor 
compared with strut II and III. Strut III has the fewest amount of unburned 
hydrogen at the outlet. The differences in loss in total pressure between the 
three struts are relatively small. 
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Fig. 5. Normalized hydrogen mass flux and total pressure over the normalized 
channel length for strut I to III, respectively 


4 Investigations of Wall Injectors 

In this section a scramjet combustion chamber with wall injection is investi¬ 
gated. Thereby the geometry and flow conditions are based on the HyShot 
flight experiments performed since 2001 by the University of Queensland to 
investigate supersonic combustion in flight. Since the beginning of the pro¬ 
gram four launches have been realized (HyShot I-IV) in which supersonic 
combustion has been achieved at HyShot II (2002) and HyShot III (2006). 
Additionally in the year 2003 and again in 2007 an identical model of the 
HyShot combustor has been investigated in the High Enthalpy Shock Tunnel 
in Gottingen (HEG) [21], Thereby several freestream conditions, equivalence 
ratios and angles of attack have been investigated. Fig. 6 demonstrates a 
sketch of the combustion chamber. The flow approaches an intake ramp (363 
x 100 mm) with a Mach-number of about 7.8. Due to the induced shocks 



Fig. 6. Sketch of the combustion chamber with wall injectors 
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300 mm 



Fig. 7. Computational domain for the numerical simulation of the combustion cham¬ 
ber with wall injectors 


at the leading edges of the intake ramp and the combustion chamber (300 x 
75 x 9.8 mm) the flow reaches a Mach-number of 2.65 in front of the wall 
injectors. To ensure a relatively homogenous flow field in the combustion 
chamber the boundary layer is bled and shocks are by-passed through the 
floor gap indicated in Fig. 6. Again the numerical simulation has been per¬ 
formed in two steps: a two-dimensional simulation for the intake ramp and 
a three-dimensional simulation for the combustion chamber and the diffuser. 
The computational domain (Fig. 7) covers a slice from the middle of a wall in¬ 
jector to the symmetry line between two portholes using symmetry boundary 
conditions Thereby the grid of the computational domain has about 1.500.000 
cell volumes. Tab. 2 summarizes the inflow conditions for the numerical sim- 


Table 2. Inflow conditions for the simulation of the combustion chamber with the 
lobed strut injector 



air 

wall injector 

pressure p [bar] 

1725 

258135 

temperature T [K] 

221 

250 

velocity u [™] 

2328 

1250 

Mach-number Ma [-] 

7.8 

1 

mass fraction Ytt 
h 2 

0 

1 

mass fraction Yo 
u 2 

0.23 

0 

mass fraction Y 

0.77 

0 

*[§] 

422.78 

4.24 
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hydrogen injection 




hydrogen injection 


Fig. 8. Calculated H2,H20 and OH distribution (from top to bottom) of the HyShot 
scramjet combustor, respectively 


illation of the combustion chamber. Fig. 8 shows the calculated H 2 , H 2 O and 
OH distribution of the HyShot combustion chamber with hydrogen injection 
into air. A lifted flame is obtained and combustion starts about 60 mm down¬ 
stream of the hydrogen injection. Fig. 8 shows in the beginning a relatively 
homogeneous hydrogen distribution in vertical direction. The majority of the 
unburned hydrogen at the end of the combustor (about 200-300 mm distance 
downstream of the combustor entrance) is located at the upper wall. The 
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Fig. 9. Normalized hydrogen mass flux and total pressure over the channel length 
for the HyShot scramjet combustor, respectively 


H 2 O-distribution demonstrates that combustion starts at the outer regions of 
the area covered by hydrogen. After ignition the flame has a circle-like shape 
and turns to a relatively homogeneous distribution at a distance of about 
200 mm from the combustor entrance. Further downstream the main combus¬ 
tion zone is located at the upper wall. Similar observation are obtained from 
the OH-distribution. Fig. 9 describes the normalized hydrogen mass flux and 
normalized total pressure along the channel length, respectively. Again the 
ignition delay can be identified by the rapid decrease of the hydrogen mass 
flux taking place at a distance of about 60 mm downstream of the fuel injec¬ 
tion. At the combustor exit nearly no unburned fuel is left. The total pressure 
shows two dips: the first due to the injection of the hydrogen at about x = 
0.04 m and the second due to the combustion process at about x = 0.11 m. 
The total pressure loss at the end of the combustion chamber is about 64% 
compared to the value at the combustor entrance. 


5 Performance 

Previous investigations of the performance of TASCOM3D on the NEC SX- 
8 [22] have shown that the performance in case of combustion simulations 
is quite good. Thereby a huge part of the whole computational time is con¬ 
sumed by a subroutine calculating PDF averaged chemical production rates 
and source term Jacobians. Due to the fact that chemistry is a local phenom¬ 
ena, a good vector lengths, quantity of MFLOPS and vector operation ratio 
are reached in the corresponding routine. Another huge part is consumed by 
the implicit solver used which is required because of the high numerical stiff¬ 
ness caused by combustion. The implicit solver is vectorized along diagonals 
(planes in 3D) to prevent data dependencies. The length of these diagonals 
depends on the size of the block and hence is dependent on the geometry 









164 M. Kindler et al. 


which could yield to relatively short vector lengths in blocks with few cells. 
In this paper additionally the MPI - parallelization of TASCOM3D has been 
reviewed and tested in order to achieve good performance in parallel runs. 
Therefore a testcase (i.e. two dimensional channel flow of an inert gas) with 
about 2.000.000 cell volumes has been calculated by a one block serial run 
and parallel runs with two, four and eight blocks, respectively. Tab. 3 shows 
the number of CPUs, the total CPU time consumed by the simulation, the 
averaged vector length and vector operation ratio as well as the MFLOPs 
reached for 1000 iterations. For the serial run TASCOM shows an excellent 


Table 3. Perfomance parameter for simulation with 1, 2, 4 and 8 CPUs 


CPUs [-] 

total CPU-time [s] 

aver, v.len [-] 

v.op. ratio [-] 

MFLOPS [-] 

1 

1008 

243.9 

99.2 

6133.8 

2 

1079 

242.7 

99.09 

5731.8 

4 

1165 

234.8 

99.01 

5307.7 

8 

1473 

205.8 

98.80 

4200.4 


performance. The averaged vector length and vector operation ratio reaches 
nearly the maximum values and the total CPU-time is 1008 s. With increasing 
number of CPUs the total CPU-time raises up to 1473 s for 8 CPUs (which is 
a factor of 1.46). The main reason can be found in the averaged vector length 
which becomes shorter with increasing number of CPUs. The total number of 
cell volumes per CPU becomes smaller due to the fragmentation of the grid. 
Hence the diagonals for the implicit solver are shorter and thus the vector 
length. Therefore the number of cell volumes per CPU should be sufficient 
high in every simulation in order to reach a good performance. The problem 
of the shortening of vector lengths would be less dramatic in three-dimensional 
test-cases (which have not been investigated in this section) because the vector 
length significantly increases while the implicit solver vectorizes along planes 
compared to diagonals. Nevertheless the performance of the parallelization 
has to be investigated more detailed and further improved. 


6 Conclusion 

The scientific inhouse code TASCOM3D has been used to investigate different 
configurations of scramjet combustors. One investigated combustor configura¬ 
tion uses lobed strut injectors while the other configuration uses wall injection. 
In case of the lobed strut injection the effect of different amounts of circula¬ 
tion on the combustion process is investigated additionally. The code has been 
tested on the NEC SX-8 and has shown to have a good vector performance. 
Also the implemented MPI - parallelization has been tested on the NEC SX-8. 
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It has been shown that the number of cell volumes per CPU should be as high 
as possible in order to reach a good performance. 
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